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pg , Abstract: We study the inflationary dynamics in a model of slow-roll inflation in warped throat. 

Inflation is realized by the motion of a D-brane along the radial direction of the throat, and at later 
stages instabilities develop in the angular directions. We closely investigate both the single field 
potential relevant for the slow-roll phase, and the full multi-field one including the angular modes 

CN ' which becomes important at later stages. We study the main features of the instability process, 

[~^. ■ discussing its possible consequences and identifying the vacua towards which the angular modes are 

^D ' driven. 
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1. Introduction and summary 

Inflation is our most widely accepted paradigm of the very early universe that enables us to explain the 
observed properties of the cosmic microwave background anisotropics. Present data are well described 
by the simplest inflationary model, consisting of a single, slowly rolling light scalar field M . However, 
data are also consistent with more complex models, for example characterized by non-canonical kinetic 
terms, or involving the dynamics of multiple scalar fields during inflation. Future observations are 
expected to enable us to distinguish among the various scenarios, provided that we can compute 
distinctive observational consequences of different setups. 

From the theoretical point of view, building models of inflation provides both opportunities and 
challenges. Inflation gives a unique opportunity to probe high energy physics since, in its most natural 
realizations, the characteristic energy scale is well beyond the current or planned particle acceler- 
ators. Also, inflation is highly sensitive to its ultraviolet completion: higher dimensional operators 
contributing to the inflationary potential play a crucial role in determining inflationary dynamics. This 



implies that inflation is able to probe properties of the theory that underlies a given model. On the 
other hand, this also means that a detailed knowledge of the setup under consideration is necessary, 
when embedding models of inflation in fundamental theories like string theory (or its supergravity 
limit). Usually, one consequence of this ultraviolet sensitivity is that a successful inflationary setup 
requires careful tunings of the available parameters to avoid the corrections due to higher dimensional 
operators, which can spoil the delicate inflationary dynamics and lead to the so-called 77 problem g, [|. 

Despite these challenges, there have been several attempts to embed consistent models of inflation 
within string theory: see Ref. B for representative examples, and Ref. B for recent reviews. The 
abundance of moduli fields in the string compactifications, both closed and open, in principle provides 
us a wide range of possible inflaton candidates. On the other hand, precisely due to the fact that many 
moduli are usually involved in a given string model, there is the danger that light moduli would interfere 
with the inflationary process. This is due to the fact that light moduli can gravitationally couple with 
the inflaton candidate (s), generically leading to the aforementioned rj problem. Fortunately, over the 
past ten years new methods have been provided to stabilize undesired moduli within string theory, 
by means of fluxes and non-perturbative effects, starting with the seminal work of Giddings, Kachru 
and Polchinski Q. However, explicitly calculable setups are generally scarce where various required 
ingredients such as the background geometric fluxes or non-perturbative superpotential 1 are known. 

Given these reasons, the warped deformed conifold [pi holds a rather special place as an ideal 
playground for string inflationary model building, leading to the framework of brane inflation in warped 
throats. To be specific, following the original work (lQ] , one considers mobile D3-brane(s) moving in 
such a background geometry, that are attracted by an anti D3-brane located at the tip of the conifold. 
The setup is sufficiently well understood that the metric, the background geometric fluxes and the 
moduli stabilization effects are known in detail, and the potential governing the entire inflationary 
trajectory can be constructed explicitly. Moreover, the parameter space of such a model is also rich 
enough that semi-realistic inflationary trajectories can be found and compared with observational data. 
Finally, D3-brane inflation in warped deformed conifold is a rare example among string inflationary 



models which has a holographic dual description M, |9|, 11 . This offers us a new perspective on the 
various contributions to the D3-brane motion from dual field theory, and allows us to employ the 
powerful and highly developed computational techniques of gauge/string duality. 

It has been shown that, by suitably tuning the ingredients responsible for stabilizing moduli, it 



is possible to find examples of inflection point inflation in this scenario |12, 13 . Inflation occurs 
in regions of the warped throat in which different competing forces act on the moving D3-brane, 
compensating each other in such a way that the resulting potential is sufficiently flat around an 
inflection point (see Refs. |1J, [lj, [16| for subsequent developments along these lines). In the original 
model, a Kuperstein-embedded D7-brane on the warped conifold ]17| was considered, and the non- 
perturbative gaugino condensation on its worldvolume contributes to stabilize all the moduli besides 
the inflaton. A sufficient number of e-foldings of single field, slow-roll inflation can then be obtained, 
which can be geometrically interpreted as a D3-brane moving along the radial direction of the warped 
throat. All the moduli that do not take part in the inflationary dynamics are made sufficiently massive 
by fluxes and non-perturbative effects. 

In this work, we revisit such an explicit framework, but instead focus on the so-called the Ouyang 
embedding |18|, p.9fl . The main motivation here is to enrich the properties of the inflationary dynamics, 
in particular obtaining a framework in which more than one field play a role in the inflationary process. 
Indeed, in Ref. |l3f (see also Refs. pO, |2l[]), it was anticipated that, for this embedding, moduli fields 

1 We however also note that, there have been interesting proposals to geometrize the non-perturbative superpotential 
in the four dimensional effective field theory, i.e. to replace them by explicit ten-dimensional geometric fluxes M pf. 



associated with the D3-brane motion along the angular directions of the warped manifold become 
tachyonic towards the tip of the throat. This property is very interesting for realizing a scenario 
in which the scalar fields associated with angular directions take part in the inflationary dynamics. 
Indeed, it is known that the evolution of the curvature perturbation after horizon exit, in multiple 
field scenarios, can have important consequences for the spectrum of the curvature perturbation and 
non-Gaussianity. The dynamics of light angular fields in brane inflation, in particular in the context 
of DBI models with non-canonical kinetic terms [£2[ |23j, have been widely studied over the past years, 
starting from Ref. |24j, mainly due to the fact that they can lead to a peculiar pattern of non-Gaussian 
spectrum of the curvature perturbation [£3| |5| . 

In this paper we mainly focus on the slow-roll inflation. We show that an expansion of more than 
60 e-foldings can be obtained in this model, by again realizing inflection point inflation. Namely, a D3- 
brane can slow-roll along the radial direction of the throat, spanning for a sufficiently long range before 
the angular directions become tachyonic. The amplitude and spectral tilt of the power spectrum of 
the curvature perturbation, produced during the epoch of slow-roll inflation, are compatible with the 
present data. The instability develops along the transverse angular directions when the inflaton reaches 
a region near the tip of the throat, and may make inflation ending a la hybrid inflation |2|, [26|. All 
the process is fully under control from the supergravity point of view. We present the necessary tools 
to follow the dynamics explicitly, and estimate the masses of the fields involved along the inflationary 
trajectory. For our particular model of inflation, we find that angular directions become very light, 
and eventually tachyonic, only towards later stages of inflation. The angular directions then roll down 
the potential towards their true minima. The fact that the dynamics remain single field for most 
of the inflationary trajectory renders the predictions of our particular model similar to the ones of 
standard slow-roll inflation. On the other hand, the mass of the angular mode during inflation is never 
significantly larger than the Hubble scale. This can render the instability process particularly slow, 
with important observational consequences for the power spectrum as discussed in other contexts in 
Ref. [£7j (see also Ref. |2^] for recent analysis of the dynamics of cosmological perturbations due to 
the waterfall field in hybrid inflation). 

To conclude this introduction, it is important to mention that hybrid inflation scenarios were 
considered in previous realizations of D-D inflationary models, as discussed in some of the papers of 
Ref. B. In those contexts, inflation ends when an instability of the tachyon field connecting the D-D 
system arises. This can lead to formation of cosmic strings at the end of inflation, which is a subject 
that has raised a large research activity starting with Ref. |2J| (see for example Ref. |3(| for a review). 
An advantage of our model is that, as we mentioned, the development of instabilities along the angular 
directions can be described within a supergravity framework, without involving string theory tachyons. 
On the other hand, as we will discuss, after the instability process takes places the system does not 
immediately fall into a global minimum. Consequently, it is not possible to make definite statements 
about the production and stability of topological defects at the end of inflation in our setup without 
relying on more detailed analysis. 

The paper is organized as follows. In Section 0, after introducing some basic facts of brane 
potentials in the warped deformed conifold, we present the explicit form of the potential in the Ouyang 
embedding. In Section pL we discuss the kinetic term of a mobile D3-brane and identify the properly 
normalized angular fields, and briefly address features of the angular masses near the tip of the throat. 
In Section El we show that, with an explicit example, that an expansion of 60 e-foldings of slow-roll 
inflation is possible in the Ouyang embedding. In Section 0, we give an analysis of the multi-field 
potential of both radial and angular fields, and provide the necessary formulae for further study of the 
multi-field dynamics. In Section El we conclude concisely discussing interesting issues to be analysed 



in the future. We relegate calculational details to three appendices. 

2. Warped conifold as an inflationary playground 

2.1 Parametrizing the potential and multiple field inflation 

In generic string compactifications, a smooth warped deformed conifold throat can easily develop in 
a bulk compact Calabi-Yau manifold near the conifold singularity. Such a deformation was done in 
the context of IIB string theory by the backreaction of localized imaginary self dual (ISD) three form 
fluxes B]. To study the D3 trajectory in such a throat, an elegant systematic way to parametrize 
the potential experienced by a D3-brane was recently given in Refs. [|8|, [ll|, which can be summarized 
succintly as 

V B3 = T 3 $_ , (2.1) 

where T% = l/[(27r) 3 (a') 2 ] is the D3-brane tension, and $_ = e A — a is a combination of the throat 



warp factor and the five form field strength given in Ref. 11 . The mode <&_ characterizes the 
perturbation away from the ISD background such as warped deformed conifold where $_ = G_ = 0. 
Here G_ = (*6 — i)G3 i s the IASD component of the complex three form flux. At low orders in 
perturbative expansion, D3-branes only couple to $_ field, and do not couple to metric or dilaton 
fluctuations, thus giving the simple form ( |2.1| ). The equation of motion of $_ can be derived explicitly 
from IIB supergravity action g as 

„8A+0 

V 2 $- = ^^|G_| 2 + 7e 4 + e - 4 - 4 |V<iM 2 +S local , (2.2) 

where e^ = g s is the string coupling, TZi is the four dimensional Ricci scalar and Si OC ai denotes the 
localized sources. The Laplacian operator V 2 is defined with respect to the unperturbed background 
metric. As the warped throat is attached to the bulk Calabi-Yau at some large radius njv, $- 
captures the bulk perturbations such as distant supersymmetry breaking fluxes or quantum effects. In 
the language of holography, at these large radius, ultraviolet perturbations are packaged into the so- 
called u non-normalizable modes" . Moreover, $_ can also receive contributions from the " normalizable 
modes'" of supergravity fields, which encode small radius, infrared perturbations. A particular example 
is the perturbation due to the D3 at the tip of warped deformed conifold which induces the D3-D3 
Coulomb attraction j31f| . In other words, $_ can in principle parametrize the potential for the entire 
trajectory of a D3-brane in the warped deformed conifold throat. 

The parametric solution to (2.2) have been extensively investigated in Refs. |8| pT[ . The results 



can be expressed in a simple form as 

Vm(<t>) = ^ D3 -D3(^) + E C *^ W0 Al , (2-3) 

i 

where c% are constants and cj> is the canonically normalized scalar field describing the radial motion of 
the D3-brane in the deformed conifold. The parameter "J = {0±, 02,4>i> < t ) 2,i>} collectively denotes the 
five different angular coordinates of the conifold 2 , and Hi(ty) is expressible in terms of the angular 
harmonic functions. We have isolated the D3-D3 Coulomb interaction as it is a normalizable pertur- 
bation, while the polynomial series consists of the contributions from non-normalizable perturbations. 
It should also be noted that there can be additional 0- independent constants added to Vd3(^>) due 



2 See Appendix W for the explicit metric and angular coordinates of warped deformed conifold. 



to the coupling with the four dimensional Ricci scalar IZ4 which can also receive contributions from 
distant sources. In Refs. || p| , much of effort was devoted to enumerating the discrete spectrum 
of the scaling dimensions {A^}. The angular harmonic function H^ty) is related to Aj and can also 



be computed in principle from the expansion of (2.2) 3 . The values of the expansion coefficients a 



however are model dependent, and can only be specified when quantities such as particular moduli 
stabilization effects, supersymmetry breaking fluxes and other microscopic quantities are known. 



The parametrization (2.3) makes it apparent that the radial and angular motions of a D3-brane 
are coupled in warped throat, and offer rich landscape for inflationary model building. On one hand, to 
obtain a single field inflationary model where <f> plays the role of the canonical inflaton, it is necessary to 
restrict to a special trajectory in the parameter space where all the angular modes are either stabilized 
at their minima or decoupled from </>: a good example was presented in Ref. Jl3[ ]. We, on the other 
hand, should expect multi-field inflation being rather generic in the warped throat. If we are to realize 
this and construct an explicit potential, it is necessary to study how D3 angular modes $ couple to 
the radial mode </>. In particular, as we demonstrate in an example in Section after taking into 
account moduli stabilization, the angular masses can change sign for certain values of </>, and become 
tachyonic. We have thus a situation where single field inflation is connected to a system with multiple 
field dynamics. 

2.2 An explicit inflaton potential: Ouyang embedding case 

We consider the motion of a D3-brane on a warped deformed conifold, whose geometrical properties 
are summarized in Appendix KL One of the basic quantities that characterize the geometry, and that 
plays an important role in what follows, is the deformation parameter e. We include a D3-brane at the 
tip of the cone, and we take into account stabilization effects that are needed for providing masses to 
undesired light moduli. In this setup, we calculate the potential experienced by the moving D3-brane. 
In this and the following sections we mostly present our main results, while details can be found in 
the appendices to which we will refer in due course. In order to calculate the D3-brane potential, we 



do not directly solve (2.2). Instead, we follow a supergravity approach, that allows us to apply well 
established techniques and results on the properties of the warped deformed conifold with embedded 
D7-branes. 

More specifically, we work in the framework of the KKLT moduli stabilization mechanism |32[ , 
where in order to stabilize the Kahler moduli 

p = a + i X , (2.4) 

whose real part corresponds to the overall volume of the compact Calabi-Yau space, we consider the 
superpotential W(z a ,p) consisting of two contributions 

W(z a ,p) = W + A{z a )e- ap . (2.5) 

Here, z 1 , • • • z 4 parameterize the complex coordinates of the internal manifold. The first term Wo = 
f G3AQ3 is the perturbative Gukov-Vafa-Witten flux superpotential Q, which, at least in principle, 
can stabilize the complex structure moduli and the dilaton-axion combination. Without loss of gen- 
erality, we shall assume Wq £ M.~ . It is well known that the tree level Kahler potential for p exhibits 



3 To be more specific, the computations in Refs. [H, hlfl were done in the asymptotic, singular conifold limit. To 
extend the analysis into the deformed conifold region, we expect logarithmic corrections to the scaling dimensions {A^}, 
and modifications to the angular harmonic function Hi(^/) as the (7(1) subgroup of the SU(2) X 5(7(2) X (7(1) singular 
conifold isometry group is broken down to discrete Z2. 



no scale structure and leaves p unfixed. A mechanism for stabilizing p is therefore to include the 
non-perturbative gaugino condensate on a stack of space-filling D7-branes (or a Euclidean D3-brane) 
wrapping a holomorphic four cycle in the Calabi-Yau space, as appears in the second term of ( |2.5| ). In 
the presence of a mobile D3-brane, the one-loop determinant A(z a ) picks up dependence on the D3- 
brane position moduli, which appear through the holomorphic D7 embedding function f(z a ) = |34[] 4 . 
We can rewrite A(z a ) as 



A(z a ) = A 



f{z a ) 

/(0) 



l/n 



(2.6) 



Here Aq is a complex constant whose exact value depends on other stabilized complex structure moduli, 
and n is the number of D7s (or n — 1 for a Euclidean D3) giving the gaugino condensate (or instanton 
correction), which also enters in the definition of the exponent a in (2.5) as a — 2n/n. In other words, 
the non-perturbative gaugino condensate on D7-branes not only stabilizes Kahler modulus p, but also 
generates potential for the mobile D3-brane, echoing an earlier result of Ref. J35|. This also fits well 
with the general parametrization ( |2.3| ), as it was also discussed in Ref. S that gaugino condensate on 
D7-branes can act as local sources for IASD flux, and make contribution to $_. 

In this paper we focus on the so-called Ouyang embedding in warped deformed conifold [R8l n9| , 



/K) =n-wi, 



(2.7) 



where w\ = {z\+iz2)/\/2 and p is a complex constant 5 . Note that \p\ heuristically measures the depth 
which D7-branes enter the conifold throat. The D7 embedding ( ^j| ) breaks the 577(2) x 517(2) x 17(1) 
isometry of the singular conifold down to 6 U(l) x £/(!), and can generate potential for the D3 angular 



modes associated with the broken isometries. The explicit potential for the w-embedding (2.7) can be 
obtained by substituting (2.5), (2.6) and (2.7), along with the warped Kahler potential Q describing 
the kinetic terms for p and D3-brane 



K 2 K,(z a ,z a ,p,p) = -3log[p + p--yk(z a ,Z a )} = -3logU(p,p,z a ,z a ) , 



(2. 



into the standard Af = 1 supergravity expression of the i^-term scalar potential. Here n 2 = M™ = 
8ttG, 7 = <7o'7 , 3/(3M|j) with <7o being the stabilized value of a at the tip of the warped deformed 
conifold, and k(z a , z a ) is the geometric Kahler potential for the deformed conifold given in ( |A.22|). The 
detailed calculations are similar to those for the Kuperstein embedding [[17] as given in Refs. |0, |l6| , 



and required inverse deformed conifold metric is given in (A.23), (A.24) and (A.25). After some 
calculations, we can obtain Vp as 



V F = AQ 



l/n 



6 _ 6ae -Mrr-V(2n) 



IA.I 



A 



•X/n-l 



V' 



coth tC 



e 2 coshr 
a^np 2 k" 



V 



(2.9) 



4 The complex coordinates defining the deformed conifold are given in Appendix pU 

5 The quantities Wi (i = 1, •••4) correspond to an alternative coordinate system for the internal manifold under 
investigation. For the relation between {ui a } and {z a }, see Appendix [A[ 

6 The residual U(l) X U(l) isometry gets further broken to U(l) in the deformed conifold. For the subtleties of 
supersymmetric D7 embeddings in warped deformed conifold, see Ref. lid]. 
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(2.10) 
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We can see that Vf depends on the volume modulus a, the radial coordinate r and various combi- 
nations of the D3-brane coordinates {wi + WT 7 \wi\ 2 ,W2 +WI, \w2\ 2 ,wiW2 +iiW}- The coordinate 
t is associated with the radial coordinate r, frequently discussed when analysing the singular warped 
throat, by r 3 = e 2 coshr: see (A. 26). In the limit A(wi) — > Ao, Vf reduces to the scalar potential 



considered in Ref. [ho|. We stress that similar computations have been done in the singular conifold 
limit of ( |2.7| ) |20| , |2l| . Here we complete the computation for the full warped deformed conifold, as it 
will be crucial for studying the angular instability at small radius in a fully under control setup. 



As the isometries of the conifold are partially broken by the D7 embedding (2.7), some of the 
angular directions can pick up effective masses through Vp, and their values are localized along the 
extremal trajectory which satisfies 



dVp 



(2.15) 



In writing the above equation, we separate the radial coordinate r from the five angular coordinates 



of the warped conifold, collectively denoted with ^. For the embedding (2.7), two such trajectories 
exist and are given by 
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(0) _ ± _L e r/2 

1 V2 ' 
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(non-delta-flat) 


(2.16) 


(0) = ± ^ e -r/2 

1 V2 ' 


»™=^" ! - 


(delta-flat) 


(2.17) 



The detailed steps are given in Appendix M. Notice that these extremal trajectories do not fix the 
overall sign of wi and W2, but only the relative sign between the two quantities. The choice of the 
overall sign has important physical consequences that we will discuss in the following sections. The 
resultant F-term scalar potential along these extremal trajectories is given by 



V F (r,a)=A(T,a)g 1 /n (r) 
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and the functions which explicitly depend on the choice of the extremal trajectories are 
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Here, we have defined two dimensionless parameters a and (3 by 
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They have the following geometrical meaning, a measures the depth which the D7-branes extend into 
the deformed conifold, while (3 is inversely proportional to the four dimensional Planck mass Mpi, and 
hence to the ultraviolet cutoff ryy of the warped throat, i.e. (3 ~ e 4 / 3 /r 2 Jv . We therefore deduce that 
r uv ZL IJ 2 i an d equivalently a > /3. In fact, when we take into account the contribution to M-p\ due 
to the compact Calabi-Yau to which the throat is connected, we generally expect a ^> (3. Also, the 
sign choice depends on the sign chosen in the extremal trajectories. Choosing the upper sign implies 



to choose the upper sign in the expressions above, and vice versa. In obtaining (2.18), we have also 



stabilized the axion \ i n l ne Kahlcr modulus p. The scalar potential now appears as a function of 
only two variables r and a, which we also necessarily need to stabilize to prevent decompactifcation. 
Before we conclude this section, let us make some remarks. First, as noted in Rcf. [}13|| , the non- 
delta-flat trajectory represents an unstable minimum. While at large radius, the angular mass matrix 
normally has positive eigenvalues, at small radius the trajectory becomes unstable and the angular 
modes become light and tachyonic. In the next sections we shall investigate in detail such instability 
and the multi-field potential. Second, we should note that in addition to the -F-term scalar potential, 
the mobile D3-brane also experiences a Coulomb attraction coming from the D3 at the tip of warped 
deformed conifold, whose form is given by 



V, 



D3-D3 



(r,cr) 



D () 
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16^Ti(A y y 



(2.28) 



where Do = 2T^a^ and Ay parameterizes the D3-D3 separation. This provides an additional contribu- 
tion to the inflationary potential. Notice that, as shown in Ref. [[10|, if the separation between brane 
and antibrane is sufficiently large, the Coulomb potential does not depend on the angular coordinates. 
Third, let us also briefly comment on the possible relevance of worldvolume fluxes on the D7-branes 
that can in principle induce a force on the moving D3-brane. As shown in Ref. [Q, this force is 
generically suppressed with respect to the other contributions to the D3 potential. Consequently, 
worldvolume flux effects are not expected to qualitatively change the analysis of the inflationary dy- 
namics that we are going to discuss. 

To conclude, let us once again discuss the difference between the supergravity approach adopted 
in this section and in the rest of the paper, and the more direct approach introduced in the previous 
section based on the dynamics of the field $_. The two approaches should provide the same results, 
as long as one neglects the effects of the compact Calabi-Yau manifold to which the throat is attached 



at large radius. The method of Section 2.1 is in principle suitable also for taking into account the 
effects of the bulk Calabi-Yau that are instead hard to describe within the supergravity approach of 
this section. On the other hand, as we will discuss in due course, the phenomena that we analyse in 
this paper are not particularly sensitive to these bulk effects, and our results are expected to remain 
correct also when we include these contributions. 

3. Angular stability analysis 



Naively, after we further stabilize the volume modulus a in (2.18) at a given value of r, we obtain a 
single field potential in r. However for this statement to hold for the entire deformed throat, we need 
to ensure that the broken angular isometries are much heavier than the radial mode and effectively 
stabilized at their minima along the entire extremal trajectories. In this section, we investigate this 
issue for all broken isometries by examining their effective masses which depend on the radial direction. 
The first step is to identify the canonically normalized fields associated with the angular coordinates. 
As we discuss in Appendix 0, when considering small displacements of the angular directions from the 
extremal trajectories, it is convenient to assemble the angular degrees of freedom into three variables 
that we call P, Q and R. We concentrate on fields representing small displacements from the extremal 
trajectory since, as we will discuss in more detail in what comes next, the angular masses during 
slow-roll inflation are large and positive. So angular modes remain localized nearby the extremal 
trajectory. Only at later stages of slow-roll inflation, as we will see, angular modes become light and 
large displacements from the extremal trajectory can occur. We will discuss the consequences of this 
fact in Section [|. 



The kinetic term of a mobile D3-brane can be deduced from the Kahler potential ( |2.8| ) as 

£kin = — ^ TT^6- (3- 1 ) 



Here, ds§ denotes the pull-back of the warped deformed metric ( A.16 ) and ct*(t) is the solution to the 
stabilization condition 

aCVp + vw) 



da 



0. (3.2) 



Note that using cr*, we can write o"o = c*(r = 0). We can obtain the analytic expression of 0v(t) for 
small r region as in the case of the Kuperstein embedding: see Appendix O. But for general r we 
have to solve the stabilization condition numerically. Considering small angular displacements around 



the extremal trajectories, the pull-back of the metric is given by ( B.65 ) as 



,4/3 



dsi\ =^rK{r) 



dr 2 + UP 2 



cosh' 



(5) 



dR A + sinlr ( - ) dQ' 



2 ~~ K "' { 3K 3 (t) 
From this form, we can canonically normalize the radial direction by defining a scalar field cj>(t) as 



(3.3) 



=( & r _*i 



(i -" y/ 6 >e ' / Vo i^roV^') 



o"o 



'22/3Vana ' P1 io #(r') V ^*(0 ' 



(3.4) 



where the function -RT(t) is given by (A. 17). Also we rescale angular fluctuations including constant 
factors as 



p =J ±T 3 e^ 3 P = 2 1 / 3 3 1 /\ ^- 
3 V anao 



R = V27^e 2/3 R = 2 1/3 3 2/ 
Q =^/2J\e 2/3 Q = 2 1/3 3 2/ 







anao 
anao 



MpiP, 

MpiR, 
Mpi Q , 



(3.5) 
(3.6) 
(3.7) 



so that they pick up mass dimension 1. We can then write the kinetic term Lagrangian as 



c^ = — 2 (d<t>) - 2 ^ (t) 



if 3 (r) 



-I- cosh 



<9i? +sinli 



9Q 



(3.8) 



where r should be regarded as a function of (^ by inverting relation (3.4), i.e. r = t(0). Kinetic terms 
with field-dependent coefficients, as the previous ones, have been used in the past in the context of 
multiple field inflation. They arise when discussing models characterized by a non-trivial metric in 
field space: see Ref. ]37| for some literature. We will use this form of kinetic terms in our numerical 
analysis in Section H. 

Let us now discuss the signs of the mass eigenvalues of the angular fields P, Q and R in the 



small t limit. A similar analysis was performed in Refs. JL3|, |20| in the singular conifold limit, and the 
results showed that some of the angular directions become generically tachyonic at small values of the 
radial direction. Here we complete the analysis to the case of warped deformed conifold, and analyse 
in detail the range of parameters leading to instabilities. 



In the limit r — > 0, the two different extremal trajectories, (2.16) and (2.17) coincide. After some 
calculations (see Appendix [j| for details), we can find that the mass eigenvalues for the angular modes 
P, Q and R, that we denote respectively with Xp, Xq and Ar, are given by 



X P = =p 2a 
Aq=0, 



Mo 



l/n— 1 r 



2s-3+— (IT a)" 2 



X R =2{lTa)X P +Aa 



zAQl' 11 - 1 ( (i 



5/3 l 



(3.9) 
(3.10) 
(3.11) 
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At this point, we can understand whether the eigenvalues are negative or positive, when approaching 
the tip of the cone (r — >• 0). Instead of doing the full systematic analysis, here we concentrate on one 
representative case that allows us to qualitatively understand under which conditions an instability 



can arise. By choosing the lower sign in Xp, with uplifting ratio s (C.4) being large enough, it is easy 
to make sure that both Xp and Xp are positive at the tip of the warped throat 7 . We focus here on 
this case, that will be the one analysed in more detail in the following sections. On the other hand, 
since Xq — > as r — > 0, suppose it approaches zero from below: since it is positive at large r, by 
continuity, this means that it changes sign in an intermediate region, and develops an instability. A 
sufficient condition to ensure that the eigenvalue approaches zero from below is to demand that it has 
negative first derivative along t, in the limit r — > 0. A simple calculation provides 



dX Q 
dr 


T=0 


2a AG 1 '"- 1 
n 


(2s- 


-3)(l + a)-a- 


4a (^ 

5/3(1 + a) \ 


a 
"5 



(3.12) 
When the ratio a//3 is sufficiently large, this is negative, but a//3 is generally much larger than one: 



see the discussion after ( 2.27 ). We see explicitly that Xq becomes negative for sufficiently small r 
while it vanishes at r = 0. These analytic considerations will be confirmed by the numerical analysis 
of the following section. 

4. Slow-roll inflation in the Ouyang embedding 

In the previous section, we found that there appears an instability at least in one or more angular 
masses along the extremal trajectories for natural choices of parameters. On the other hand, by 
suitably tuning the parameters, it might be possible to realize slow-roll inflation that lasts for more 
than 60 e-folds well before the instability develops. In this section, we show that this is indeed the 



case. We specifically focus on the lower branch solution of the non-delta- flat direction (|2.16| ), i.e. 
Wi = — ee r / 2 /\/2 and W2 — te~ T I" 1 j \[2. The discussion on the application of the delta-flat trajectory 
to the multi-field inflation is postponed in the next section. The choice of the parameters is summarized 
in Table [FJ 



a (3 ns \W \ A 

0.00444 0.00157 8 1.1254 4.64 x lO^Af^ O.OO&M^ 

Table 1: The parameters used in the numerical computations. 

The single field potential along the radial direction r is shown in Fig. |l| in terms of the canonically 
normalized field 4>(t). Here we plot the potential in terms of </>(t)/(^, where 



2 — u II «JI/f2 



:-! 



4/3 



.^- Wi- (4-1) 

(ftp corresponds to the maximal value for <p in the warped throat. We see that there is a point where 
the second derivative of the potential with respect to vanishes. Slow-roll inflation can be realized 
around this inflection point. In order to understand whether inflation can be realized or not, we define 



7 Notice that this does not imply that Xp and Xp are always positive at r > 0. 
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the slow-roll parameters 8 



H 
IP 



V 



M, 



pr 



v 



(4.2) 
(4.3) 



These slow-roll parameters are shown in Fig. [2|. As expected, e remains small but \rj\ quickly becomes 
large once the field moves away from the inflection point. 
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Figure 1: The full potential versus normalized canonical field 4>/4>in which runs from to 1. We set Mp\ = f . 
In the left panel, we simply show the entire potential, while in the right panel we concentrate on the inflecting 
region, where after the instability of Q develops we show the potential with a dotted line. 
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Figure 2: The "potential" slow-roll parameters (left) e and (right) r\. Note that while \r\\ ~ 1 when the 
instability of Q develops, e remains smaller than 1: at the moment of instability, we have e ~ 3.58931 x 10~ 4 . 

Next we consider the masses of the angular directions. In Fig. y, we show Xp, Xq and Xr. For our 
choice of parameters, we find that all the angular directions become unstable at small radius. However 
we should note that once instability develops, we cannot trust anymore the analysis of the mass matrix 



8 Notice that we write the first slow-roll parameter as e to distinguish it from the deformation parameter e used 
elsewhere. 
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around the extremal trajectory in Appendix |B| as the angular displacements need not to be small since 
angular fields start to roll towards their true minimum. Thus we only study slow-roll inflation before 
one of the angular directions becomes unstable: in this case the first among the angular directions 
whose mass squared vanishes is Q. From this point on, we need to consider multi-field dynamics, 
which is the subject of the next section. 



X 





0.00 0.05 0.10 0.15 0.20 0.25 0.30 

<PI<t>n 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 



Figure 3: The angular masses (left panel) Xp, (right panel) (dotted line) Xq and (solid line) Xr. The 
difference between Xq and Xr is only noticeable near the tip region. 



We can estimate the number of e-folding number from this instability point as 

2a 



N 



3/3 



(4.4) 



2ancTQ 

which is shown in Fig. H As we can see, the number of e-folds blows up near the inflection point. That 
is, we expect an indefinitely large number of e-folds if the classical initial condition is such that the field 
starts from very near the inflection point. Of course, in reality quantum fluctuations will push the field 
away from the inflection point. In any case, we learn that it is not difficult to realize slow-roll inflation 
that lasts for more than 60 e-folds before instability develops along one of the angular directions. An 
interesting, separated issue is the problem of initial conditions that can lead to inflation: we will not 
discuss it in this context, but see for example Ref. [ p8| . 

Also, we calculate the power spectrum of the curvature perturbation V-n and the corresponding 
spectral index n-jz, under the slow- roll approximation given by 

^=24^' (4 - 5) 

n-jz =1 - 6e + 2j] . (4.6) 

In Fig. o, we show log 10 V-r, and n-ji versus the number of e-folds. For our choice of parameters, we 
obtain V n = 2.44096 x 10" 9 and n K = 0.936381 at N = 60, which are well within the 2a range of the 
current observations ]l|. 

The observational predictions from slow-roll inflation is very similar to those in the Kuperstein 
embedding m& . This suggests that the inflection point inflation is quite generally realized regardless 
of the choice of the embeddings. In fact, as discussed in Section 0, it was suggested that the D3-brane 
potential can be expressed as 
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Figure 4: The number of e-folds N versus cj>/cj>^ counted from the instability point. As we approach the 
inflection point N diverges, which indicates that we have indefinitely large N near the inflection point. 
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Figure 5: (Left panel) log 10 Viz and (right panel) n-n versus N. At N = 60, the values of Viz and n-n arc 
close to the observed values. 





where \I/ is the angular coordinates, c$ are constants and A, are given by 

A, 



1, 2 , 2, 2 , 3, 



(4.8) 



In our setup, we find that the potential along the extremal trajectory is well fitted by the polynomial 
form 



V 



V Q 



1 + 0.090 



( ^ 



VMpi 



0.108 



Mpi 



3/2 , , 

- 2.637 f 



VM P1 



•6.316 



Mpi 



5/2 



3.678 



(4.9) 
where Vq = 2.250 x 10 -16 . This constant term arises from the uplifting contribution that is necessary 
to obtain a quasi de Sitter solution. This parametrization emphasizes the fact that inflation is obtained 
around an inflection point region, that results from a delicate cancelation among the various terms 
in the previous expansion. Notice that although the previous expression provides a potential for the 
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canonically normalized inflaton along the radial direction <j>, the dependence on the angular coordinates 
is contained in the coefficients Hi(ty). Their form can be deduced from the explicit scalar potential 
derived earlier, which also gives mass eigenvalues Xp, Xq and Xr along the extremal trajectories. 
The knowledge of these quantities is necessary to understand whether or not the angular masses are 
heavier than the Hubble parameter during inflation: if not, angular directions are not stabilized on 
their extremal value and their dynamics can have interesting observational consequences. 

To estimate the effective masses for the angular directions, as can be read from (jO| ), we notice the 
kinetic terms induce a mixing between the radial field 4> and the angular fields P, Q and R. Precisely 
speaking, we need to properly take into account this non-trivial field space metric to analyze the masses 
of the angular fields and we will explain the procedure in more detail in the next section. However, 
during inflation, the motion along the radial direction is suppressed by the slow-roll parameter. Thus 
at the leading order in the slow-roll approximation, we can assume that the radial field is constant, as 
well as the coefficients of the kinetic terms. We are then interested in the effective masses, that take 
into account the non-trivial r-dependent coefficients in the kinetic term, as 



> _ X P aU ,...', T-l, 



^^^V^^' (4 ' l0) 

m Q,eff M 2 { 22/3 3 4/ 3/3 K{t) sinh 2 (r/2) > ^^> 

m 1 = *« ™ lAll (4 12) 

m fl,eff M ^ 22/3 3 4/ 3/3 X(T ) cosh 2 (r/2) ' V' l£ > 

For a given 0, these masses characterize the behaviour of the isocurvature perturbations along the 
angular directions. We found at N = 60, their ratios to the Hubble parameter are 

"=5.31529, (4.13) 



m- 

P.cff 



H 2 



<2 



^£-=2.00828, (4.14) 



o 
R,e« 



,.,, =1.99503, (4.15) 

H z 

and the behaviours of these ratios are shown in Fig. g. We learn that the effective masses of the 

angular fields are comparable to the Hubble parameter during inflation: they are stabilized on their 

extremal values until the instability region is reached, at which the fields become tachyonic. 

This conclusion is not unexpected. In supergravity, the typical mass of the moduli fields during 

inflation is of 0{H). This is the reason of the r\ problem that we mentioned in Section ttc only by 

carefully tuning the parameters, one can obtain a flat potential for the inflaton field in some region 

of parameter space. In our framework, we can also understand more explicitly why the masses of 

the angular modes turn out to acquire these values during inflation, using the large r expressions for 

XpQp. These quantities, in this limit, are found to be 

X P ^±^(l T ae^) - 2(1/ "~ 1) ae ^ Uaa + - - 3«e°« (l T ae^) ^ 
n \ J y 2 \Aq\ \ ) 

±2ae T/2 (l T ae 



1 =f ae T/2 
3V3 n jj 



r .I) ^ 3 2/3^ e r/3 



(4.16) 



X Q ,X R ^ (l T ae^ 2 ) X P + ^-^e^ {l T ae^y 2 ^ lM f-l + ^^A . (4.17) 
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Figure 6: The plots of (left panel) m~ ft/H 2 an d (right panel) (solid line) m| rc/H 2 and (dotted line) 
m~ a/H 2 versus 4>/4>ij.- Note that we assume H is completely dominated by the potential, which should be a 
good enough approximation near the flat region of the potential. 



With our choice of parameters, the stabilisation condition for a can be approximated by (CIO), which 
gives <t* ~ (To even for large r. Then substituting ( C.10| ) into the above expressions, we find that the 
term proportional to aa cancels out where aa* ^$> 1 in order for the supergravity approximation to be 



valid. Comparing these expressions with the potential (2.18), we can estimate Xp t Q t p as 



A-P,Q,R = V, 

ana* 



(4.18) 



where Mp^q^r are functions of ae T ^ 2 , a 2 e T / 3 //3 and s. In the example we studied here, ae T / 2 <C 
1, a 2 e r / 3 //3 -C 1 and Np,q,r — 0{ae T / 2 ) at N = 60. Then the effective masses can be estimated as 



m 



P.Q.R.cS 

H 2 



O 



(3 



(4.19) 



This is generally larger than one as a > /J, but they are not significantly large with our choices of 
parameters: the mass eigenvalues during inflation turn out to be of the order of the Hubble scale. 
This can have important consequences for the field dynamics right at the end of inflation. Indeed, in 
different contexts, mainly motivated by hybrid inflation, it has been shown that when the waterfall 
fields have not too large masses the instability process might take longer than one Hubble time to 
complete |27|. This implies that slow-roll motion along the radial direction can continue, while the 
angular modes roll, not necessarily too rapidly, towards their true minima. It would be interesting 
to investigate these issues in the present context, although the analysis has to take into account the 
specific form of both the potential and the non-canonical kinetic terms. 

Let us finish this section discussing whether the results of this section can be modified after taking 
into account the effects of the bulk Calabi-Yau attached to the throat at large radius. As shown 
in Ref. p9], the main point to notice is that bulk contributions to the angular potential give rise 
to suppressed by powers of the warp factor. Then, modifications of the angular potential due to 
bulk effects are generally subleading with respect to the potential produced by moduli stabilization 
effects on the throat. On the other hand, bulk effects could play a role in the infrared region where 
the angular instabilities develop, and where the overall force on the angular modes due to moduli 
stabilization effects vanish. But bulk contributions are subdominant away from this particular region, 
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being suppressed by the warp factor. Consequently, we expect that their net effect is at most to induce 
a tiny shift of the position at which the instability occurs along the radial direction. 

5. Explicit multi-field potential from warped conifold 

In our previous analysis, we mostly focused on studying the homogeneous dynamics of the inflaton field 
along the extremal trajectories. On the other hand, since towards the later stages of slow-roll inflation 
an instability develops in the angular directions, in this region the masses of the corresponding fields 
become small, comparable to that of the radial direction, and eventually tachyonic. To follow the 
dynamics from this stage onwards, we need to consider in more detail the properties of the angular 
fields, significantly deviating from their extremal values. With this purpose, in this section we study 
more closely the potential governing the angular fields that develop an instability towards the end 
of inflation, as well as the corresponding kinetic terms. This analysis represents the necessary first 
step for discussing in detail the field evolution in this system at the homogenous level as well as the 
dynamics of fluctuations. Simple considerations, based on the properties of the potential, allow us to 
obtain important information about the instability process. In particular, we show that the instability 
connects one of the extremal trajectories followed along during inflation to another, e.g. non-delta-flat 



trajectory to delta-flat or non-delta-flat one. This is studied in Section 5.2 where we also discuss 
the fact that the actual dynamics of the system is characterized by non-standard kinetic terms which 
contain cross terms involving derivatives of the fields. 

5.1 Alternative parametrization of the angular fields 

The fields P, Q and R that we used in Sections g and || are given by combinations of the angular 
coordinates a*, ft and ji in the warped conifold, using the so-called Euler-Rodriguez parametrization 
discussed in Appendix y| Recall that, in Section^, we found that the angular field P remains heavier 
than Q and R during the most interesting phase of inflation, and does not take part in the instability 
process 9 . Thus we set it at its extremal value P = in the following discussion. From its definition 
(B.33|), this gives a,\ + a 2 = 0. We will restrict ourselves to the case ol\ — a 2 — 0, which greatly 



simplifies the following calculations. We then concentrate on the remaining angular modes Q and R in 
the previous discussion. We have seen in the previous section that the masses of these modes become 
tachyonic at almost the same point along the inflationary region, see Fig. [3| For this reason we are 
interested in determining the complete potential that describes these fields together, as well as the 
radial direction. As can be read from (B.31) and (B.32|), they are actually formed by a combinations 



of two pairs of angular variables. It turns out that, while for small perturbations around the extremal 
trajectories the combinations of /3i, ft>, 7i and 72 giving Q and R are the most convenient to deal 
with, this is no longer true when describing arbitrary displacements from the extremal points. To 
describe this last case, it is useful to define another basis of fields. 



From ( |B.5| ) and ( B.12 ), we can easily find that starting from a point along the extremal trajectories, 



u»i and u>2 are given by 



(0) I, || (0) . 1 1 . 1 1 / 71 • Pi \ ( 72 . Pi , , , ... 

w\ =w\ ; cos|xi|cos|x2| -w\ ; sin I xi I sin 1*2 1 I - — --%- — -I I - — - -%- — - ) , (5.1) 

w 2 =u4 0) cos|xi|cos|x 2 | +u4 0) sin|xi|sin|x2| (r^T+^r^T ) ( T^l + i T^l ) > ( 5 - 2 ) 

VlXil \Xi\J \\X2\ |X2 1 



9 Notice that this fact holds only for the specific inflationary model we have considered. There could be other 
inflationary trajectories, with different choices of the parameters, for which the behavior of P is different. 
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with Xi given by ( B.ll ) as X f = fit + if ■ The actual range of the quantities ft and 7^ depends on the 



extremal trajectory one considers, as discussed in Appendix B.3 (see also the discussion in the next 
section), but we can take them positive. Along with \i an d X 2 ; it is natural to dehne a field £ as 



a + 6=2£. 



-1 { 7i 
£1 = cos - — - 

Xi 



£ 2 = COS 1 



72 
Ixal 



A 
Ixil 

' 1X2 1 



(5.3) 
(5.4) 

(5.5) 



It is convenient to consider £ rather than £1 and £2 separately since, as can be read from ( [5.1| ) and 
(p^), only the combination £1 + £ 2 appears in the potential. With these new fields, and using (2.16) 



and (2.17), we can write W\ and w-i as 



± —=e T/2 L_±i—=e r/2 sin | Xl | sin | X2 | sin(2£) , (non-delta-flat) 

«'. = <> V e (5-6) 

± —=e- T,2 Lt ± i^=e T/2 sin | Xl | sin | X2 | sin(2£) , (delta-flat) 

v 2 v 2 



U> 2 



T -pe t/2 LJ; ± z— e r/2 sin | X i| sin | X2 | sin(2£) , (non-delta-flat) 

V 2 v2 

T ^=e T/2 L: ± i-^-=e' T/2 sin | Xl | sin | X2 | sin(2£) , (delta-flat) 



V2 



V2 



(5.7) 



where we have defined 



L% =cos| X i|cos| X2 | +exp(±r)sin| Xl |sin| X2 |cos(2£), 
L ± =cos|xi|cos| X2 | -exp(±r)sin| X i|sin| X2 |cos(2£). 



(5.8) 
(5.9) 



In the limit of vanishing X i, X 2 and £, we recover (2.16) and (2.17). Note that these expressions have 
an interesting property: starting from one extremal trajectory, we can move to another by continuously 
varying X i, X2 and £: for example, as will be demonstrated momentarily, from the lower branch of the 
non-delta-flat extremal trajectory obtained by X i = X2 = £ = 0, we can obtain both branches of the 
delta-flat one by taking | X i| = | X2 | = n/2 and £ = or £ = w/2. We will discuss some consequences 
of this fact in what follows. 



5.2 Potential and kinetic terms 



Having found w\ and u> 2 as (5.6) and (5.7), we can straightforwardly calculate the scalar potential. 



After some computations, we can find that the form (2.9) does not change, but the functions Q, C 
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and T> become different. The results are 



Q = 



(lTae T / 2 Llf + 1 



[sin|xi| sin|x2| sin(2^)] , (non-delta-flat) 



{l^ae- T / 2 L±) +a 2 e r [sin|xi|sin|x 2 |sin(2^)] 2 , (delta-flat) 



(5.10) 



2 tanhT I ±ae r / 2 L7 - a 2 e T LZL\ 



+ a 2 e T [sin |xi| sin |x2 1 sin 



(2£)] 2 } , (non-delta-flat) 



2tanhr[=Fae _r/2 L+ + u 2 er T L+_L% - a 2 e T [sin |xi| sin |x 2 | sin(2£)] 2 j , (delta-flat) 



(5.11) 



V 




1 — sechr 



x < tanhr + 2csch(2r) 



[sin|xi|sin|x 2 |sin(2£)] 




coth ', 



(cos|xi|cos|x2|) 2 + (sin|xi|sin|x 2 |) 2 + 7:COshrsin(2|xi|)sin(2|x2|)cos(2^) 



/ e _2 e , 2 2 

-sechr -=-£_ H — ;j— £_ + coshr [sin |xi| sin |x2| sin(2£)] 



(non-delta-flat) 



1 — sechr 



P -L-J 



e 



[sin|xi|sin|x 2 |sin(2£)] 




cothr 



x < tanhr + 2csch(2r) 



(cos |xi| cos |x 2 |) 2 + (sin |xi| sin |x 2 |) 2 + -coshr sin(2|xi|) sin(2|x 2 |) cos(2£) 



/ e ' , 2 e' _ 2 2 \ 

-sechr ——LI + — L_ + coshr [sin |xi| sin |x2 1 sin(2^)] I > . (delta-flat) 



(5.12) 



Plugging these expressions into (2J)), we can obtain the full potential including angular directions. 
This is the complete form of the potential that is needed in order to study the dynamics of the angular 
and radial directions in the instability region. While in Section we focussed on the lower branch of 
the non-delta flat trajectory, the previous formulae allow us to also obtain the potential in all other 
cases. In the limit of small Xij or equivalently of small Q and R, this potential coincides with the one 
we used in Appendix |B| to determine the eigenvalues Xq and Xr. 

The potential is periodic along the angular directions xi, X2 and £, as it is clear given our expres- 
sions for to i and W2- There are "extremal points" , defined as the positions for which derivatives of the 
potential along \i an d £, vanish. These extremal points coincide with the extremal points along the 
directions Q and R: namely they provide the same non-delta-flat and delta-flat trajectories determined 
in Appendix H. They are shown in Table 0, but others that are periodically identifiable with these 
ones are also extremal points. 



IXil 

\X2\ 



NDF, L NDF, U NDF, U DF, U DF, L 






■n 


tt/2 


tt/2 


IT 





tt/2 


tt/2 


tt/2 


tt/2 





tt/2 



Table 2: The values of xii X2 an d £ for extremal points. "NDF" and "DF" stand for non-delta-flat and 
delta-flat, and "U" and "L" for upper and lower branches respectively. For example, the first column indicates 
the lower branch in the non-delta flat trajectory, the case on which we focussed in Section 4. 
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These observations are sufficient to extract important qualitative information about the outcome of 
the instability process. For definiteness, let us consider the explicit example we studied in the previous 
section, which corresponds to the first column of Table |2j. For the choice of parameters summarised in 
Table HI, all of the extremal points have at least one unstable angular mode except for the upper branch 
of the delta-flat direction (|xi| = \\2\ = t/2), which corresponds to the true vacuum of the D3-brane 
along the angular directions. Until the point at which instability develops, all the angular modes are 
at the origin. Nearby the instability region, the non-delta-flat trajectory becomes a local maximum 
(or saddle point) for the angular modes. The angular masses become light and then tachyonic, and 
the fields \i and £ start to roll down the potential. In this case, the range of the variables are found 
to be £ € [0, 4tt] and \% G [ — 7r /2, tt/2] (or equivalently \\i\ € [0, 7r/2]), as discussed in the last part 



of Appendix B.3. After a possibly complicated dynamics, they reach another extremal point of lower 



energy, at which all the squares of the mass eigenvalues are positive with our choice of parameters. 
This corresponds to the upper branch of the delta-flat trajectory, the fourth column of Table 0. The 
actual result depends on the detailed dynamics of the system. On the other hand, our knowledge of 
the potential is sufficient to show that the angular field dynamics, starting at the point of instability, 
connects the initial non-delta-flat with one of the delta-flat trajectories. As a representative example, 
we plot in Fig. p] the potential depending on both radial and angular directions, taking for definiteness 
Xi = X2 and keeping £ = 0. The plot clearly exhibits an instability along the non-delta-flat direction 
(xi = X2 = 0), and we can more easily see that the potential connects the false vacuum corresponding 
to the lower branch of the non-delta-flat trajectory with the true vacuum corresponding to the upper 
branch of the delta-flat trajectory (|xi| = |X2 1 = 7r /2) in the angular directions. 
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Figure 7: The multi-field potential as a function of the radial direction and a particular choice of angular 
coordinates \i — X2 and £ = 0. We can see that other stable trajectories exist at |xi| = [xa] = 7r /2 next 
to the unstable one at \i = X2 — 0. The potential is periodic along the direction xi, with period n. This 
is consistent with our observation that \xi\ £ [0,7r/2]: see the last part of Appendix B.3.We used the other 
parameters the same as in Table 1. 

As we have seen from the previous arguments, the angular instability connects different extremal 
trajectories in the warped deformed conifold. Depending on initial conditions, after the field falls from 
one trajectory to the other, a non-trivial dynamics continues along the radial and, possibly, along the 
angular directions. This implies that, at the later stages of slow-roll inflation, the inflaton field does 
not generically fall into a global minimum, but into other extremal points that further evolve due 
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to non-trivial motion in the radial direction. In this setup, consequently, it is not possible to obtain 
a quantitative understanding of the process of formation and evolution of topological defects at the 
end of inflation, based on the symmetries of the potential. Instead, a detailed numerical analysis is 
necessary, which takes into account of the post-inflationary evolution along the radial direction. 

The analysis of the actual dynamics that develops at the instability is further complicated by the 
non-standard form of the kinetic terms characterizing the fields. In order to obtain them, one can 
follow an identical procedure to the one discussed in Section @ and Appendix pL After straightforward 



calculations, we can find that the form of the kinetic term is given by (3.1) where the pull-back ds^ 



given by (A. 16), with 

g™ = sin(2 Xl ) dfr + cos(2£) sin(2 X2 ) dfc + 2 sin(2£) d X2 , (5.13) 

g^ =2d X i + sin(2£) sin(2 X2 ) dfc - 2 cos(2£) d X2 , (5.14) 

g^ = sin(2 Xl ) dfr - cos(2£) sin(2 X2 ) c% 2 - 2 sin(2£) d X 2 , (5.15) 

g<-*> =2 d Xl - sin(20 sin(2 X2 ) d& + 2 cos(2^) d X2 , (5.16) 

g^ = [1 - cos(2 Xl )] dfr + [1 - cos(2 X2 )] d& . (5.17) 

The kinetic terms are then characterized by the cross terms involving derivatives of the fields. These 
cross terms have very small coefficients when considering small deviations around the extremal tra- 
jectory. But they play an important role when large angular displacements are considered, which is 
expected during the instability process. The range of the angular fields, in the previous equations, 
depends on the initial trajectory one is considering, as discussed in the first part of this section. 

In this section, we have provided all the necessary tools that are needed in order to study the 
dynamics of the instability process at later stages of inflation. We have obtained important qualitative 
information about the field dynamics and the instability process. These ingredients are the starting 
point to perform a more detailed analysis of the field evolution after the instability develops. For 
example, it would be interesting to study how the annihilation of D3-D3 branes occurs, the role of 
the angular directions in this process, and whether moduli stabilization induced by the presence of 
D7-branes can have some effects in this phenomenon. Indeed, in Ref. |39| it has been shown that near 
the tip of the throat the force between brane and antibrane can acquire significant angular dependence. 
Another important issue is to investigate possible formation and evolution of topological defects after 
inflation ends. We hope to return to discuss these issues elsewhere. 

6. Future Directions 

In this work, we have considered a model of D3-brane inflation in a warped throat. Different forces 
act on the D3-brane, due to the presence of an anti D3-brane at the tip of the throat, and of a 
Ouyang-embedded D7-brane in the warped deformed conifold. We have shown that these forces 
compensate each other, and can obtain 60 e-foldings of slow-roll inflation around an inflection point of 
the potential along the radial direction. At later stages of single field slow-roll inflation, the angular 
directions develop instabilities after which more complicated dynamics follows. Our work provides a 
new explicit example for obtaining slow-roll inflation in warped throats, other than using Kuperstein 
embedded D7-brane |l3J, |lq| . Moreover, it is an example where angular directions play an active role 
for determining the inflationary dynamics while the model is fully under control from the supergravity 
point of view and the complete potential for the relevant modes is known. 

There are various topics, motivated by our study, that deserve further analysis. We have shown 
that the masses of the angular modes remain comparable to the Hubble scale during most of the 
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inflationary period, becoming light towards the later stages of slow-roll inflation. The fact that the 
angular masses are not much larger than the Hubble scale can lead to a situation in which the slow- 
roll motion continues along the radial direction for few e-foldings, even after the development of the 
instability in the angular directions. When the slow-roll conditions are violated during the instability 
process, however, the DBI nature of the D-brane action shows up, since the brane velocity increases. 
This fact can also play an important role in determining the detailed dynamics governing the instability. 
It would be interesting to study the consequences of these features for the spectrum of the curvature 
perturbation, taking into account the presence of the non-canonical kinetic terms for the fields involved, 
and the particular form of the potential. 

Another important direction is to study the formation and stability of topological defects produced 
after inflation. As we have argued, the instability process connects different extremal trajectories. 
Consequently, after the onset of instability, the system does not reach a global minimum, and further 
motion of the brane along the radial direction takes place. The study of implications of this brane 
motion, for a possible development of topological defects, needs a detailed numerical investigation. 
Finally, it would also be interesting to understand whether the energy stored in the inflaton field can 
be transmitted into the matter sector, while the angular fields roll off their false minima. We hope to 
return to analyse these issues in the near future. 
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Appendices 

A. Properties of warped deformed conifold 

A.l Parametrization of the angular coordinates 

Here we collect some facts concerning the various coordinates parameterizing the deformed conifold 



following Ref. 40 . It can be defined as a submanifold C 4 via the complex embedding equation 



e^ q : 



2 = e 2 . 



(A.l) 



where {z*} are complex coordinates in C . They can be expressed in terms of a 2 x 2 matrix as 



W 



—w 3 w 2 \ 


1 


( z 3 + iz 4 


z 1 — iz 2 


—w 1 w 4 J 


~T2 


\z 4 +iz 2 


- z 3 + iz 4 



(A.2) 



The deformed conifold relation (A.l) can then be written as 



A ,2 ,3 ,4 



det W = w 1 w 2 - w 3 w 4 ■■ 

The complex embedding coordinates of deformed conifold { 

of the real global coordinates {r £ R , ip £ [0, A-k] , 9\^ £ [0, w] , 01,2 £ [0, 2tt}}, S = r + iip as 

, / - i , ■■ i + "1 \ I 9\ + 02 \ • • , / - \ I "' "" r ',' | I "I + ":■' 

cosh I — | cos I | cos I | + 1 smh | — | cos | | sm 



(A.3) 



z 2 =e 



z 3 =e 



z 4 =e 



— cosh — cos 
2 



i sinh ( — ) cos 



cos 



2 

'1 + 02 



— cosh I — ) sin I — ) cos I + i sinh 



. , ■ • , . | h + #2 \ . / 01 - 02 

C0SM 2" J Sln I — 2 — ] Sm l — 



, z } can also be expressed in terms 

(A.4) 

(A.5) 
(A.6) 
(A.7) 



1 — <72 \ . / 91 — 92 

2 y sin I — ^ — J sin I — ^ — 



..._.. . ji - C7 2 \ 91-92 

i sinn — sm cos 



One can readily verify that the constraints (A.l) and (A.3) are satisfied. We can also encode the 

(A.8) 



coordinates (A.4)-(A.7) compactly in the matrix form as 

W = L c W s Rt , 
where the 2x2 matrices W s , L c and R c are given by 



/ 



W s 



^l 2 \ 

V2 



Vvf 



-r/2 



L c 



cos I — exp 



sm | — I exp 



/ 

kV'i + 01 



2^i- 0i ) 



sm I — I exp 



cos I — I exp 



-Ml ~ 0l) 



-2^1+01) 



(A.9) 



(A.10) 



J / 
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and similarly for R c with the subscript 1 replaced by 2. L c and R c are two independent SU(2) group 
elements, which generate the entire deformed conifold by acting on the special point W s . Notice that 
the coordinates if)± and if>2 only appear in the combination if) = ipi + rp2- 

A. 2 The metric of the warped deformed conifold 

It is convenient to define the following basis of one forms |j, J41J] , 

sm8id(j)i — (cos if) sin $%d4>i — sva.if)d6%) 



7C 1 ) 







V2 




(2) 


d9i 


— (sin if> sin 6*2^2 + c °s if}d9 2 ) 








V2 




(3) 


— wa.6\d(f>i + (cosip sin 62(^)2 ~ sin 


^<Wa) 






V2 




(4) 


_d$i 


+ (sin if) sin 9 2 d(f>2 + cos if>d92 ) 





<] 



V2 
( 5 > =dip + cos 9\d(f>\ + cos 6*2^02 • 



(A.H) 

(A.12) 

(A.13) 

(A.14) 
(A.15) 



The explicit metric of the six dimensional internal part of the geometry, the deformed conifold, can 
be expressed as 



4/3 

d4 =^rK(r) 



1 



3K 3 (t) 
+ sinlr 



dr 2 



(i) 



,(i) 



(5) 






9 



(3) 



( 9 -) ; 



where 



ff(r) 



[sinh(2r) - 2r 



1/3 



2 1 /3 s inhr 
We can also write the full ten dimensional warped metric as 



dsi 



2A ^r 1fil/ dx^dx v + e- 2A ^ds\ 



(A.16) 



(A.17) 



(A.18) 



The warping is obtained by turning M units of F3 flux through the A-cycle of the deformed 
conifold, and — K units of H3 flux through the dual -B-cycle, when we put a ultraviolet cutoff of the 
warped throat at e~ ~ 1. The warp factor is given by the expression [B| 

e -4A(r) = 2 2/3( 5sMa ')2 e -8/3 /(T) ^ 



where 



7(r) 



dx XC0thx - l ^u^ Q^l/3 



sinh a; 



[sinh(2a;) - 2a;] 



(A.19) 



(A.20) 



The internal, deformed conifold metric ( A.16 ) can be obtained from the so-called "little" Kahler 
potential k(z a ,z^) as 

g a p = d a dpk. (A.21) 

As the angular directions of the warped deformed conifold arc isomctries, they do not appear explicitly 
in the little Kahler potential. It thus only depends on the radial coordinate r. An expression for k(r) 
is H 

F 4/3 

(A.22) 



k (r)^ e ¥N J o T dr'[smh(2r')-2r'} 1/3 , 
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where without loss of generality we set the integration constant to zero. For the calculation of the 
F-term scalar potential, the unwarped inversed Kahler metric in terms of the holomorphic embedding 
coordinates {w a } is also needed and is given by 



k^ 



k" 



5/3 _K&B a' P 



R 



If* 



i? a/3 + cothr (^-- cothr ) I" 



aB 



e 2 ct,w a wP + ci,w a wP' 



W a W@ + C® cl w a 'w^ 



(A.23) 

(A.24) 
(A.25) 



Here c^, is a 4 x 4 matrix whose only non-zero entries are c\ — c\ — — 1 and c| = c| = 1. Notice that 
the metric ( A.25 ) is indeed six dimensional as the deformed conifold constraint w\W2 — W3W4 = — e 2 /2 



allows us to substitute away e.g. W4 and W4. 

To recover the familiar singular conifold limit, we can define the radial coordinate r as 



(A.26) 



r 3 = ^|z Q | 2 = e 2 coshr. 

a=l 

In the asymptotic limit e T — > 00, the metric ( |A.16[ ) reduces to usual conical form 

3 

ds 6 -> - (dr 2 +r 2 ds 2 T i.i) , 

where the metric ds 2 ^!^ for the five dimensional base T 1,1 can be easily deduced from ( A.16Q . Fur- 
thermore the warp factor e~ 4A ( A.19| ) acquires an AdS form 



(A.27) 



-4A(r) _ 37rg s JV(a / ) 2 
~ ^4 



(A.28) 



where N is the quantized five form flux. 



B. Extremal trajectory and angular mass matrix 

In Section g, we derived the complete potential controlling the dynamics of the D3-brane along the 
throat. It explicitly depends on the D3-brane coordinates, in particular on the angular directions. In 
this section, we determine extremal trajectories of the D3-brane that extremize the potential along 
the angular directions. These extremal trajectories, along which the brane moves only in the radial 
direction, correspond to stable points in the moduli space of angular directions, at least for sufficiently 
large radius r 3 = e 2 coshr. For small values of r, instead, extremal trajectories become unstable for 
some of the angular directions. The consequences of this fact for inflation is discussed in Sections [| 
and||. 

We denote the five angular coordinates of deformed conifold by {^i}- We then look for extremal 

trajectories that satisfy 

dV 

(B.l) 



d^ % 



0. 



The chain rule ensures that finding points satisfying (B.l) is equivalent to identifying points such that 
d(wi +wT) _ <9|ui!| 2 _ d(w 2 +Wi) _ <9|u>2| 2 _ d(w 1 w 2 + wiw 2 ) 



d^,. 



dV, 



a*,- 



C^ ? ; 



a*, ; 



0. 



(B.2) 



25 



This is because, as we can read from ( |2.9| ), the i^-term scalar potential Vf only depends on these 
combinations of coordinates of the deformed conifold. Thus, our next task is to relate the brane 
coordinates w a to the angular coordinates $ \. 

To parametrize the angular coordinate *5>i, we can rewrite the matrix W as jl3| 



Here Wo is given by 



W = LeWoRl 



,(°) „„(°) 



W 



(°) ..,(«) 



— Wi w 



(B.3) 



(B.4) 



and denotes a fiducial point, which we consider to be angularly stable. The matrices L e and R e are 
two independent SU(2) group elements. Using the generators of the SU(2), i.e. the Pauli matrices 



o"i 



(To = 



V3 



L e and R e can then be identified as 



1 

1 

-i 

i 

1 
-1 



JTi 



Re =e lT2 , 
where T± and Ti are parametrized by two real three-vectors \u%, fy, 7i} (i = 1, 2) as 

Ti = 13,(71 + 7*02 + a t a 3 = 
Denoting 



O; 



Pi + iji 
2 2 i o2 i 2 



ft - «7i 



(B.5) 
(B.6) 

(B.7) 



(B.8) 
(B.9) 



(B.10) 



(B.ll) 



a simple calculation provides 



' cos xi +«"i - 



L« = 



IXi 
sin |xi 



J (ft - *7i) 



SID 



in|xi| \ 

Xi 
sin |xi 



(B.12) 



\ i(ft+i7i)— ; — |— cos|xi| -iai — ■ — -I 
\ IXil Ixil / 



and similarly for R e . 

The two SU(2) group elements L e and R e take a point Wq on deformed conifold to the entire 
manifold. Notice that here we have adopted different parametrization of SU(2) group element (Euler- 



Rodriguez) from the one used in (A. 8) (Cayley-Kline), as the former parametrization is particularly 
suitable to study the properties of extremal trajectories. The parameters {c^ ,ft ,7,} (i = 1,2) of 
( B.12| ) and {Oi , fc , ifii} (i — 1,2) of (A. 10) are usually related non-linearly. We will discuss the 
relation between these two different parametrizations in Appendix [ 
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B.l Linear expansion: Identifying the extremal trajectory 



As we want angular stability near w a satisfying (BJ2), expanding (|B.8| ) and (BJ)) and only keeping 
the linear order, then we find 



W = (1 + iTi) W (1 - iT 2 ) = W + i (IiW - W T 2 ) 
The changes in w\ and w 2 are then given by 



(B.13) 



8w\ 



,(°) 



i(«i + a 2 )u4 + (i/3i - 71)11$ + (ij3 2 - 72) 



(0) 



(0) 

J 4 i 



<Sw 2 =i(a\ + a 2 )w 2 + (ifii + ji)w\ + (ifi 2 + j 2 )w { 3 



(B.14) 

(B.15) 



Then, it is trivial to find, by demanding the linear variations of the combinations of w± and w 2 which 
appear in the potential vanish, that iw^ ,w 2 > € K and wjj — w\ — 0. Thus, the conifold 
constraint equation now gives 



(0) (o) e 

w\ 'w 2 = -— . 



Combining with ( A.2q), which now reads 



(B.16) 



..(«) 



„(°) 



KO'+KO 5 



= e coshi 



(B.17) 



we determine the following extremal trajectories 



w 



(0) 



V2 



r/2 



u; 



(0) 



T ;r 



-r/2 



(B.18) 



and 



„(°) - 






r/2 



(0) 



T— e T/2 



(B.19) 



In the large radius limit e r ' 2 — > 00, the trajectories (B.18) and ( |B.19[ ) correspond respectively to the 
so-called non-delta-flat and delta-flat trajectories determined in Ref. [[13). While we present an explicit 
realization of slow-roll inflation using the non-dclta-flat trajectory, we will also discuss the role of the 
delta-flat trajectory for the multi-field inflationary models. 



B.2 Quadratic expansion: Mass matrix 



We now develop the tools necessary to investigate the stability of the extremal trajectories (B.18) and 
( B.19| ) along the angular directions. To find the mass matrix associated with the angular coordinates, 
we expand (B.8) and (B.9) up to second order in perturbations. Explicitly, we have 



IT = j 1 + ^T 1 - ^tA W q (l - iT 2 - \sTl 



--W + i{T 1 Wo-W T 2 ) + 



TiW T 2 - i (TfWo + W Ti) 



(B.20) 
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Using w- 



(o) - ..,(<>) - 



= along the extremal trajectories, we can explicitly calculate the combinations 



{wi + W\, |wi| 2 , w 2 + w 2 , \w 2 \ 2 , W1W2 + w 1 ii! 2 }, and find 

Wl + W =2w { ° } - [(a, + a 2 f + ft 2 + ft 2 + 7l 2 + 7 |] wf - 2 (Aft - 7l72 ) wf , 
w 2 + wt =2^ 0) - [(ai + a 2 ) 2 + ft 2 + 0* + 7l 2 + 7 |] u;< 0) - 2 (ftft - 7 i 72 ) w { ° } , 

l- 2 | 2 = (4 0) ) 2 - (^i + # + 7? + 7l) (- 2 0) ) 2 - 2 (A/3, - 7l72 ) 2 w^w^ , 
Wl u; 2 + i« =2u)[ 0) ^ 0) - 2 (ft 2 + ft 2 + 7l 2 + 7 2 ) wfM 0) - 2 (ftft - 7l72 ) "...(°) 



(B.21) 
(B.22) 

(B.23) 

(B.24) 
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As required, the departures from the fiducial point vanishes at linear order. We can thus write the 
mass matrix as 



d 2 v 



d^id^ii 



dV 



d 2 {w l +W[) dV d 2 \ Wl Y 



d{w 1 +W[) d^id^j 
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d(w2 + W2) d^id^j d\w 2 \ 2 d^id^j d{w 1 w 2 +wTm) 



d^id^, 




(B.26) 



where the linear partial derivatives of the brane coordinates do not appear since we have chosen the 
extremal fiducial trajectory. We can see that in the present angular coordinates {ai,ft, 7 i}, the mass 
matrix is not diagonal because of the terms which mix different coordinates, such as (ai + a 2 ) 2 , ftft 
and 7i 7 2- Let us rewrite the angular coordinates ft and 7$ in terms another set of real coordinates g, 
and r, as 



Then we immediately find 



01,2 = 

71,2 = 



gl±ri 

V2 ' 

r2 ±g 2 

V2 ' 



Pi + Pi + 7i + 7 2 =q\ + <& + r\ + ' 



ftft - 7172 
We can further define shorthand notations 



ql + ql 



q\ 



q\ =Q 2 , 
rl=R\ 



B.27) 

B.28) 
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B.30) 



B.31) 
B.32) 



and similarly 



a\ + a 2 = P , 



(B.33) 
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since only these combinations appears. Using these new coordinates {P, Q,R}, we obtain 
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These new coordinates therefore make the mass matrix diagonal, with rows and the columns corre- 
sponding to {P ,qi ,q 2 ,ri , r 2 }, as 
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where 



X P 

X, 



,(°) 



dV 



d(wi + wi) 

■ Q = -2(«,i°)+«4 0) ) 



«,\ 



(«) 



<9U 



9(^2 + W2) 



9U 



9(wi + wi) 



dV 



+ 



dV 



d\wi\ 



2 [w\ 



(0) 



If 



2[u4°>-u;f 



+ e z 



dV 



9(w 2 + tl) 2 ) 
dV 



d\w 2 \ 



-2 to. 



+ 2e 2 (1 -coshr) 



(0) 



(B.40) 
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(B.42) 



It is consequently straightforward to substitute the derivatives of the potential, as well as the 
expressions for w\ in the previous formulae. Explicitly, we can calculate the angular masses associated 
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with them as 

=F2a (e T/2 V Wl+m - e~ T/2 V W2+W2 -)\ , (non-delta-fiat) 



X P = 



Xn=< 



=F2a (e t/2 V Wi+ — - e T / 2 V W2+m )\ , (delta-fiat) 

X P + 2a 2 (2V WlW2+W[W2 - — e T V\ Wl \2 — e~ T V\ W2 \2) (non-delta-fiat) 

+ [2a 2 (V\ Wl \2 +V\ W2 \2) - 4a 2 coshrK, lU , 2+ww =F 2a (e r/2 V W2+W2 - - e- T / 2 V Wl+m: )]\ 

X P + 2a 2 (2V WlW2+wrm2 - - er T V\ Wl \2 - e T V\ W2 \2) (delta-flat) 

+ [2a 2 (Vj m |2 + V\ W2 \2) - 4a 2 coshrK, lU , 2+w1E j =F 2a (e~ T / 2 T4, 2+ ^ - e r/2 V w - L+ - m )] | Q 



(B.43) 



(B.44) 



X R = l 



X P + 2a 2 {2V WlW2+wpuH - e T V\ Wl \2 - e~ T V\ W2 \2) (non-delta-flat) 
- [2a 2 (V\ Wl \2 + V\ W2 \2) - 4a 2 cosh tV Wi W2 +wTm T 2a (e T/2 V W2+W2 - - e~ T / 2 V Wl+m )] \ Q 
X P + 2a 2 (2V WlW2+w[W2 - - er T V\ Wl \2 - e T V\ W2 \2) (delta-flat) 
k - [2a 2 (V\ Wl \2 + V\ W2 \2) - 4a 2 coshrK, lU , 2+wlI ^ =F 2a (e~ T / 2 K, 2+ ^ - e T/2 V Wl+m )] \ Q 



(B.45) 



where 



n 



■A^l/n-l 



B - 3ae a °^g- 1/{2n) - ?— - (K coth rC + 1l V a ) - K cothr 
\Ao\ ny 



v - A c x l n - X - 

y w 2 +w 2 — vq • , 

n smn r 



Vi 



«'i 



v, 



W-L+Wi VQ 



A l/ n -l 

— C/n cothr 



JC + TZ X ( --cothr 



Vi 



\w 2 \ 



V, 



WlW2TWl^2 



n 

— Q a cschr 
n 



and 



I + cothr I — — cothr 



JC + TZ X I — -cothr 



_ k' _ 3 sinh(2r) - 2r 



W -T77 = 



k" 4 sinh 2 r 



n ( T ) = 



anp?^k' 



ftx(r) 



1 



= 2a 2 cosh tK x (t), 
3 2 / 3 1 [sinh(2r) - 2r] 2 / 3 



" awyk" 2 4 / 3 /3 



sinh 



(B.46) 

(B.47) 

(B.48) 
(B.49) 
(B.50) 

(B.51) 
(B.52) 

(B.53) 



These results are the basic tools for studying the stability of inflationary trajectories along the angular 
directions. 

B.3 Euler-Rodriguez versus Cayley-Kline parametrization of SU(2) group 

We now discuss in more detail the connection between the Euler-Rodriguez and the Cayley-Kline 
parametrizations of SU(2) group elements. In particular, we are interested in the relation between 
the matrices L e and R e , and the matrices L c and R c . The analysis here also leads to the derivation 
of the kinetic terms used in the multi-field analysis in the main text. 
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Recall that W is obtained separately by using two different parametrizations, (A.8) and (B.3) 



Since, the expression for Wq (B^) depends on which extremal trajectory one chooses, also the relation 
between L e and L c depends on this choice. From the similarity transformation between Wq and W s , 
we can find that one possible realization of the relation between L c and L e is 



L c = < 



L e ia2 , (non-delta-flat, lower branch) 

L e ia\ , (non-delta-flat, upper branch) 

L e 1 , (delta-flat, lower branch) 

L e d3 , (delta-flat, upper branch) 



(B.54) 



and precisely the same for R c and R e . But other realizations work as well. For the remaining part 
of this appendix, we focus on the lower branch of the non-delta flat trajectory: this is the case we 
explicitly considered in Section 0, and the other cases can be treated in an identical way. The relation 
between the angles 9i, ipt and <f>i used in the Cayley- Kline representation and a*, /3i and 7^ in the 



Euler-Rodriguez representation is obtained from the first equation in (B.54), and reads, both for i = 1 
and 2, 



cos 



cos — sin 



sin — cos 



sm 



(to 


+ & 




2 




2 

-4>i 




2 

- 4>i 



= ~7» 



\Xi\ 
sin \Xi 



\Xi\ 
cos|xi|, 

sin|x»| 
\Xi\ ' 



(B.55) 
(B.56) 
(B.57) 
(B.58) 



By means of these relations, we see that the extremal trajectory corresponding to c^ = /3j = 7^ = 



is obtained, for example, by choosing 9 
around this point, we find 



(0) 



71", tp. 



(0) 



,(0) 



27T. By expanding the above expressions 



Sai = - (6<f>i - 5ipi) , 



8 fy = — - sin 



Sji 



4 0) + 4>. 



4>: 



(0) 



(0)' 



k(f>r 



(B.59) 
(B.60) 

(B.61) 



These expansions allow to obtain the relation between {P,Q,R} and {6i,<f>i,ipi}. Choosing for defi- 



nitcness ip. 



(0) 



27T, 



do) 



0, the identifications QB.31| ), p.32| ) and ( |B.33|) yield 

Q = - ^8(61 - o 2 ) , 



R 



2^/2 
1 

2V2 



8(6! + 2 ) 



(B.62) 
(B.63) 
(B.64) 
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It is easy to check that, near this extremal trajectory, expanding up to quadratic order in angular 



fluctuations, the pull back- metric (A. IS) reduces to 



,4/3 



d ~4\ = ^rK{r) 



dr 2 + 4dP 2 
3K 3 (t) 



mslr ( -) dR 2 + sinh 2 (-)dQ 2 



(B.65) 



This result plays an important role in Section 0. We have checked that the same expression holds for 
the expansion around all the extremal trajectories. 

To conclude, it is interesting to consider in more detail the case a\ = a-i = 0, that plays an 
important role in the discussion of Section p. In this case, in order to satisfy (B.58), we can choose 
ipi = 2tt + (pi. Then ( B.55| ), (B.56) and (B.57) simplify, and provide 



\X 



2 2 ' 

Pi =IXi|sin-0i, 
It =\Xt\cosipi. 



(B.66) 

(B.67) 

(B.68) 



They are valid for arbitrary values of the angles, inside their interval ranges. Since 9i e [0, w], we then 
learn that \xi\ G [0,7r/2], or equivalently xt € [ — 7r /2, t/2]. Moreover, comparing (B.67) and (B.68) 



with ( |5.4| ) and (5.5), we can make the identification £ = ip in this specific case. This implies that £ lies 
in the range [0, 4tt]. The potential is periodic in these angular variables, with periods set by the ranges 
we have just found. However, let us emphasize once again that the identifications ( B.66J ), ( |B.67 ) and 
( B.68J ) and the relative field range, are only valid starting from the lower branch of the non-delta- flat 
trajectory. Making other choices in ( B.54 ) lead to different identifications. On the other hand, it is 
simple to work out them following the same procedure as above. 



C. Stabilized value of a 

In this Appendix, we analyse the stabilisation of a in the limit r — > 0. In the absence of uplifting, the 
stabilized value of a at r = is given by the condition 



dV F (0,a) 



c)a 



= 0. 



From (2.18), one easily finds 



\ A o\ 



Then, we add to Vf an uplifting term of the form 



V D (r,cr) 



D(r) 
U 2 (r,a) 



from which we define the uplifting ratio 



_ D(0)/U 2 (0,a F ) 



1 + aap 

2 + aap 



J-(O). 



\V F (0,a F )\ 
In order to avoid a runaway decompactification, s should be of order one, 1 <s < 3 p 



(C.l) 



(C.2) 



(C.3) 



(C.4) 
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The stabilized value of the Kahler modulus after uplifting, which we call oo = a F + 5a J13J, is 
found from the condition 



where 



Then in the limit 5a / ' a F <C 1, we obtain 

5a 



dV 
~da~ 


o-o 


da 2 


A J. dvD 

do- 

of 


dV D 


2V D (^ Ja\ 




da 


CO * 


'V'**,)- 



a F 



3+[oHd*V F /da2)/2Vn}\ aF 
In the limit of large aa F , we can show that 

d 2 v F 



da 2 



2a'\V F (0,a F )\ 



Plugging this expression into flC.7j ) and using {CA), we obtain 

S 



aSa 



aa F 



(C.5) 
(C.6) 

(C.7) 
(C.8) 
(C.9) 



Then it is then easy to see that 5a <C 1 induces the following contribution to the stabilized modulus, 



3ae Q 



|Woj 



Go 



l/(2r. 



'(0) 



3 + 2aa + + acrF jr( ) + 2s . 
2 + aa F 



(CIO) 



The stabilized value of a, in the limit of r — > 0, is given by the solution of ( p. 10 ). 

The dependence of a on r is more difficult to determine, for non-zero r, and a numerical analysis 
is needed. See however Refs. |1S, n& for some analytical results. 
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